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We study the closeness of an experimental unitary bosonic network with only partially indistin¬ 
guishable bosons in an arbitrary mixed input state, in particular an experimental realization of the 
Boson Sampling, to the ideal bosonic network, where the measure of closeness of two networks is the 
trace distance between the output probability distributions. An upper bound on the trace distance 
to the ideal bosonic network is proven and also a bound on the difference between probabilities of 
an output configuration. Moreover, the upper bound on the trace distance is tight, provided that 
a physically transparent distinguishability conjecture is true. For a small distinguishability error it 
is shown that a realistic device with N bosons is at a constant trace distance to the ideal Boson 
Sampling under the 0(A'^“^)-scaling of the mismatch of internal states of bosons. 

PACS numbers: 42.50.St, 03.67.Ac, 42.50.Ar 


I. INTRODUCTION 

The Boson Sampling (BS) of S. Aaronson and 
A. Arkhipov [ij can potentially become the first oper¬ 
ating device with postclassical computational power, if 
scaled up to few dozens of bosons. It can be imple¬ 
mented with passive linear optical devices, single-photon 
sources, and non-adaptive detection of photons, cleaBy 
not enough for the universal quantum computation 
1^. A complex interference of bosonic path amplitudes 
in a unitary linear network elevates the output proba¬ 
bility distribution of the BS to the #P class of compu¬ 
tational complexity asymptotically inaccessible for 

a classical computer even to verify the result. The key 
point in optical implementation of the BS is the quan¬ 
tum indistinguishability of single photons, resulting in 
the Hong-Ou-Mandel dip 0 (see also Ref. Q). Sev¬ 
eral labs have performed the proof of principle experi¬ 
ments on small networks with few single photons [l^ 
El. Various error models have been analyzed [TsI - llSj . 
It is now clear that an experimental device with post- 
classical computational power is feasible, if setup errors 
are properly tuned [T^ . An experimental implementa¬ 
tion can now be chosen between various setups [20l - l^ . 
Though no unconditional verification method of an oper¬ 
ating device is known, there are results in this direction 
[ 2 I with an experimental confirmation [23, El . More¬ 
over, there is a conditional certification protocol based 
on many-body interference El . Recently, it was shown 
that the Franck-Condon transition probabilities in molec¬ 
ular vibronic spectra can be simulated on the BS with a 
modified input El - 

In the present work we consider how to measure the 
closeness between a realistic BS device, implemented 
with only partially indistinguishable bosons, and the 
ideal BS. We use the trace distance between the output 
probability distributions of a realistic device and the ideal 
BS as the measure of their closeness and study how to es¬ 
timate the trace distance using the information derivable 
from experiment. We give an upper bound on the trace 


distance, and show that under a highly plausible phys¬ 
ically clear conjecture the derived bound is tight. We 
employ a recently developed partial indistinguishability 
theory ( 2 ^ in our derivation of the upper bound on the 
trace distance. In fact, we derive the main results in a 
more general form, so that they can be applied to mea¬ 
sure the closeness of an arbitrary experimental bosonic 
network, with an arbitrary (mixed) input state of bosons, 
to the ideal bosonic one with the same distribution of 
particles over the input modes. Thus one can use the 
upper bound derived in the present work to quantify the 
“bosonness” of an experimental network. Moreover, the 
upper bound on the trace distance between a realistic 
and the ideal bosonic networks involves the probability 
of the maximally bunched output (when all particles are 
in the same output mode), which can be estimated from 
the experimental data. 

Here it should be noted that the purity measure of 
partial indistinguishability, proposed in Ref. [2^ - which 
quantifies the quality of a quantum interference in an 
unitary linear network (for instance, it reduces to the 
Mandel’s parameter E^ for two bosons) is not suitable, 
in general, for attesting how close is a realistic network 
to the ideal bosonic network (i.e., with the output ampli¬ 
tudes proportional to the matrix permanents of subma¬ 
trices of the network matrix). This is due to the fact that 
both species of identical particles, bosons and fermions, 
can behave in an interference experiment as the other 
species, if the particles are in an appropriate (entangled) 
input state (one can even simulate the BS with fermions) 
[SOj. For instance, the output probability distribution 
of a bosonic network with bosons maximizing the purity 
measure of Ref. E^ (thus being completely indistinguish¬ 
able particles) can be identically equivalent to that of the 
ideal fermionic one, if the input state is anti-symmetric 
with respect to permutation of internal states of parti¬ 
cles. In the latter case the bosonic network is still a per¬ 
fectly quantum network (that is why it has the maximal 
purity), but equivalent to the ideal fermionic network. 
Only in the case when the input state of bosons is just 
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a small perturbation of a state symmetric with respect 
to internal states of particles, the purity measure can be 
used to quantify how close such a network is to the ideal 
bosonic one. In general case, to measure the closeness 
to the ideal bosonic network a new specific measure is 
required, which must reflect how the internal state of 
bosons (the projection of the input state of particles on 
the Hilbert space of the degrees of freedom not affected 
by a network) is close to a nearest symmetric state, since 
only such a symmetric input state results in the ideal 
bosonic behavior [^ . 

Since there is a continuous family of input states for 
every probability distribution realizable at output of an 
unitary network [^ . it is more convenient to work with 
the partial indistinguishability matrix (J-matrix) of Refs. 
[iM , which is the unique image of all input states of 
bosons resulting in the same output probability distribu¬ 
tion in a unitary network. A linear space Vn of dimen¬ 
sion N\ for N particles can be associated with J-matrices, 
where vectors and matrices are labeled by elements (per¬ 
mutations) of the symmetric group Sn- A crucial point 
is that a boson counting measurement at network out¬ 
put is formally a measurement also in the auxiliary lin¬ 
ear space Vn- Thanks to this fact, by application of 
the standard methods we derive a uniform (over all net¬ 
works) upper bound on the trace distance of an output 
probability distribution to that of the ideal bosonic net¬ 
work, which supersedes the probabilistic bound of Ref. 
M- Under a physically transparent distinguishability 
conjecture, there is also a lower bound on the trace dis¬ 
tance, which has the same scaling in the total number of 
bosons as the upper bound, thus making the latter tight. 
The distinguishability conjecture is a restricted from of a 
highly plausible property that different J-matrices result 
in distinguishable probability distributions, when the lat¬ 
ter are compared over all possible networks. 

In the case of bosons in pure spectral states, our upper 
bound on the trace distance to the ideal bosonic network 
is proportional to a recently proposed permanental mea¬ 
sure [^. However, our bound is tighter by 1/N\ (the 
proportionality factor) for N bosons. 

The rest of the text is organized as follows. In an aux¬ 
iliary section III Al we give a brief account of the results 
of Ref. [ 2 ^ used below. In section III HI we prove Theo¬ 
rem 1 (the main result), giving the upper bound on the 
trace distance to the output probability distribution of 
the ideal bosonic network. Some mathematical details 
of the derivation are placed in appendix In section 
III Cl we consider a special case of bosons in pure inter¬ 
nal states. In section III PI we derive a bound on the 
difference between probabilities of an output configura¬ 
tion, Theorem 2 (with some mathematical details placed 
in appendix [g. In section m we study the limit of a 
small distinguishability error and give a scaling law of 
the mismatch of internal states of bosons with the to¬ 
tal number of particles for a constant trace distance to 
the ideal bosonic network. In section m we consider a 
lower bound on the trace distance and argue that under 


a highly plausible distinguishability conjecture our upper 
bound is tight. In section|V]we give concluding remarks. 


II. MAIN RESULTS 

A. Output probability formula for a general case of 
partially indistinguishable bosons 

Let us first recall the necessary details of the partial in¬ 
distinguishability theory in a unitary linear bosonic net¬ 
work (see for more details Ref. [1^). 

Consider a linear unitary M-mode network with N 
identical bosons at its input, such that each input mode 
k = receives a certain number of particles 

Uk > 0. The single-particle degrees of freedom are di¬ 
vided into two parts. One part is composed of modes 
operated on by a linear unitary network. The other 
part consists of internal states unaffected by the network 
(whose Hilbert space will be denoted by "H). A network 
is given by a unitary matrix U, which leaves the Hilbert 
space of internal states invariant, i.e., acting as follows 

M 

= ( 1 ) 

where o-lj/bl.j creates a boson in an input/output 
mode k and an internal state \(j)j) (the set of states 
(l^i),..., Ifij ),...} is a basis of "H). A general input state 
of configuration n = (ni,..., um) (we will use such vec¬ 
tor notations below) reads p{n) = 
with Qi >0, and 


N 


= 


V Tin) 




fee ,1, 


JO), 


( 2 ) 


OL — 1 


where pin) = 0^1 ki,...,kN are (generally re¬ 

peated) input modes, and |0) is the vacuum state. 

Due to permutation symmetry of boson creation op¬ 
erators, the expansion coefficients can be chosen to 

satisfy permutation symmetry with respect to a Young 
subgroup of the symmetric group Sn, consisting of per¬ 
mutations TT of internal states of bosons in each input 
mode between themselves, tt G S^ 0 ... ® = 5^, 

where Sn is the symmetric group of n objects. Such co¬ 
efficients are normalized by X)? = 1- 

For a given unitary network, probability of detection 
of an output configuration rh = (toi, ..., mM), |hr| = N 
with |to| = TOi -I-... -I- iriM, depends on the degree of par¬ 
tial indistinguishably of bosons, defined by their internal 
states. Moreover, experimental particle detection is im¬ 
perfect, prone to particle losses and dark counts, besides 
introducing an effective filtering of the internal states due 
to non-ideal sensitivity of detectors. We assume that such 
effects are small and consider the postselected case of all 
bosons being detected at the output of a network. There 
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still remains to consider the effect of non-ideal sensitiv¬ 
ity of detectors (e.g., of limited spectral width in the case 
of photons). Here we adopt the simplest model of iden¬ 
tical detectors. For identical detectors the input state 
alone determines the bound on the postselected proba¬ 
bility distribution at output of a unitary network (since 
in this case the same J-matrix, defined below, is assigned 
to any output configuration of bosons). 

The particle counting measurement is therefore de¬ 
scribed by the following positive operator-valued mea¬ 
sure {n(m), \m\ = N}, where each element reads (see 
the derivation for photons in appendix A of Ref. [T^ l 


n(m) = 


^(m 


E 


N 


Ilrj 

.a=l 


N 


N 


Q!=l 


n 




( 3 ) 

here I = {li,... ,lpf) is a set of output modes, gener¬ 
ally repeated, corresponding to mode occupations m, 
0 < Fj < 1 is sensitivity of a detector to an internal state 
\(j)j). In the case of non-ideal photon counting detectors, 
the summation over discrete internal states is replaced 
by an integral over spectral states [2^. Probability of 
detection of an output configuration rh reads 


p{rh\n) = Tr{n(TO)p(n)}. 


( 4 ) 


Inserting Eqs. and ([3]) into Eq. (jd]) we obtain output 
probability as a quadratic form, which can be expressed 
as a double sum over Sn [H, 


p{rn\n) = 


1 


N 


p{m)p{n) 


^ ^ <Ari,(T2 




criCa): 


Uk 




( 5 ) 

The partial indistinguishability matrix J indexed by two 
permutations (Ti,cr 2 G Sn (which has only up to A^! dif¬ 
ferent elements for identical detectors) is given by a trace 


= Tr 


( 6 ) 


where f = Tj\(f)j){(f)j\^ is a representation of tr G 

(Sat in , defined as 


Pcr\(t)n) ®---® 

(7) 

and is an internal state of bosons correspond¬ 
ing to an input state of Eq. i.e., pb"*) = 

N 

( 8 ) 

3 

The ideal bosonic behavior in a unitary linear network 
occurs for = 1 for R is known that this 

case realizes if and only if the internal state of bosons 
satisfies SnP^^^^'^Sn = [2^, where Sn = 


is the projector on the symmetric subspace of TL®^. In 
this case Eq. ([5]) contains the absolute value squared of 
the matrix permanent of a matrix U[n\m], built from a 
network matrix U by taking rows with multiplicities ft 
and columns with multiplicities fn. For ideal detectors 
Tij = 1 and 


p{m)pL{n) 


(9) 


In particular, the BS setup corresponds to an input with 
rifc < 1 in the dilute limit M ^ N’^. In this case probabil¬ 
ity of a bunched output (i.e., with some nik > 2) vanishes 
as 0{N'^/M) and in the asymptotic limit N ^ oo the 
output probability distribution cannot be simulated on a 
classical computer 0 , 10 . Recently it was argued that 
for M > N'^ the many-body correlations of indistinguish¬ 
able particles may not vanish in the thermodynamic limit 
[ 3 ^ . On the other hand, it is also known that for M <C iV 
the output probability distribution can be approximated 
to a small multiplicative error by a semiclassical approach 


For non-ideal detectors, Fj ^ 1, output probabilities 
p(m|n) of Eq. ([5|) do not sum to 1. To compare with 
the ideal network Eq. dH), we consider the postselected 
probabilities, when N input particles are detected at a 
network output. The postselected output probabilities 
are given by division of those of Eq. (jd]) by the detection 
probability pd = We have (see also Ref. [ 2 ^) 


^ (f , (10) 


where tr(A) = Xo- is the trace in the linear space 
Vn- 

Finally, there is a continuous family of internal states 
of bosons resulting in the same probability distribution 
at output of a unitary linear network [^ . However, it 
is highly plausible that no two different J-matrices give 
coinciding output probability distributions over all uni¬ 
tary networks 37|. Our main results, formulated below, 
do not change if the above property is not true in general 
(see also section m- 


B. Upper bound on the trace distance to the ideal 
bosonic network 

We will work in the linear space Vat, where J-matrix 
lives (for convenience, Dirac’s notations will be used for 
vectors in Vn, e.g., a vector of the natural basis, in which 
J-matrix is given by Eq. (jb)), will be denoted by |tT), 
a G Sn)- We consider how close is the postselected out¬ 
put probability distribution of a realistic bosonic network 
to the ideal bosonic network. The postselected output 
probability distribution will be denoted by p, whereas a 
particular output probability by p{rn) (dropping the in¬ 
put configuration argument ft). The output probability 
distribution of the ideal bosonic network is denoted by 








4 


a particular output probability (to) is given by 
Eq. dni). Our goal is to estimate the trace distance 

’ p) = ^ E (^) - p(^) I • (11) 

m 

The trace distance is a 1-norm based measure of closeness 
between two probability distributions (all norm based 
measures in a finite-dimensional linear space are equiva¬ 
lent). It is used in discussion of the computational com¬ 
plexity of a non-ideal BS (for details, see Ref. i)- The 
following result is found. 

Theorem 1.- The trace distance Eq. (HU satisfies 

Zl(p(®''),p)<l-p„ (12) 


with Pa defined as follows 

Pa = —Tr{f®^5wp(“*^S'jv} = —Tr{f®^p(®”*)S'Ar} 

(13) 

(the second form is due to mutual commutation oi Sn 
and and that Sat is a projector). 

Note that Pa is the postselected (for non-ideal de¬ 
tectors) probability of a symmetric internal state of 
bosons at a network input. For a symmetric input state 
we get Ps = 1. In section HTCl below it 
is shown that the only parameter in the upper bound, Ps, 
can be estimated from the experimental data, i.e., even 
when no knowledge about the input state of particles is 
available. 

Proof of Theorem 1.- The bound of Eqs. (fT^ - (fT^ eas¬ 
ily follows from the following observation. Boson count¬ 
ing measurement at a network output can he thought of as 
a “measurement” also in the auxiliary linear space Vn- 
Indeed, let us consider the trace-normalized J-matrix de¬ 
fined as follows (see Eq. (ITUl) ') 


Ja 


<J<7i ,(72 

N\pd 


^_Tr{f®^p(“‘)p i}. (14) 


(see appendix 1X1 for more details). The projector V acts 
as the identity operator on any j7-matrix; VJ = JV = 
J, due to the symmetry 

for all TT from the Young subgroup Sfi = 8^ ® ® 

Sum j as discussed in section llT^ The postselected output 
probability assumes a “measurement” form in the linear 
space Vtv 


p(to)= ^ {Zj:\j\z^), ( 18 ) 

m(Z )—'m 

where the summation is over all output mode indices I = 
(li,..., In) corresponding to an output configuration to 
( in total N\/p{m) terms). 

The rest is an elementary algebra. The independence 
of the sum ffa,cri from a (recall that depends 

on the relative permutation only) means that such 

a J -matrix has an eigenstate | s), whose expansion in the 
standard basis reads ((t|s) = l/vWI, for any permutation 
a. The corresponding eigenvalue is precisely the postse¬ 
lected probability of a symmetric input state Pa- Indeed, 
from Eq. (1141) by comparing with Eq. (USD we get 

(s|j7|s) = ^ ) i7o-i,(T2 = ^ ) Juj = Ps: (19) 

< 7 i,(T 2 <7 

where I is the identity permutation. Therefore, j7-matrix 
as a convex combination of two matrices 

J = Pa\s){s\ + {l-pa)J^^\ JW|S)=0 . (20) 

Using Eq. O to rewrite the trace distance m as fol¬ 
lows Zl(p(®'^),p) = (1 — pa)D{p^'^‘^\p*--^'>) (with cor¬ 
responding to noting that the trace distance is 

bounded by 1, we get the upper bound of Eq. HU- 

Q.E.D. 

Finally, it is notable that the upper bound Eq. m co¬ 
incides with the trace distance between a given jfT-matrix 
and the pure ^/(®‘^)-matrix of the ideal case = |s)(s|. 
We have 


This J"-matrix has all the properties of a density ma¬ 
trix in the linear space Vn: whereas the following vectors 
\Zf) G Vn (here I = {li,..., In) the output mode indices), 
in the standard basis defined by 


T>(|s)(s|, = itr|||s)(s| - } = 1-Ps 


where we have used Eq. (Eni- 


( 21 ) 


{<^\Zr) 


1 


N 




( 15 ) 


C. Upper bound on the trace-distance for pure 
internal states of bosons 


define a positive operator-valued measure in Vn- 

'^\Zr){Zr\ = r, (16) 

r 


with the projector V given as 


Pu.a' = 


fi{n] 




tr(iP) 


TV! 

/i(n) 


It turns out that, for bosons in pure internal states, the 
probability Pa defined in Eq. m is intimately related 
to the permanental measure of Ref. [3l|. Indeed, let N 
bosons have pure internal states \pa): a = 1,... ,N. We 
have from Eq. ([6]) 

N N 

J<T1,<T2 = (‘T’o-da) |r|</^(T2(a)) = (‘Pa|r|V’cr2(Ty®(a))’ 

CK—1 Q;=l 

( 22 ) 


(17) 
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thus 

N 

Pd = W_{^a\T\^a) (23) 

a—1 

and 

p-=E n !».,(.,)=( 24 ) 

a q ;=1 

where the Gram matrix built from the internal states 
reads Gap = {ipa\^pp) with |^„) = \ipa) / {{Pa\^Pa))"^ ■ 
However, our bound given by Eq. m is tighter by 
1 /N\ than the previous bound, A^! — per (G), proposed in 
Ref. im. 


D. Bound on the variation of an output probability 
from that of the ideal network 


Let us recall that the case of classically indistinguish¬ 
able particles corresponds to J" = (the maximally 

mixed j/-matrix as allowed by the symmetry = 
JV = J, see also Ref. [l^). The classical probabil¬ 
ity of an output configuration m reads 




per(|[/p[fi|m]) 

pirn) 


(25) 


where |C/p[ri|rn] is a matrix constructed from elements 
by taking rows with multiplicities ft and columns 
with multiplicities rh. We are interested in how close is 
the postselected probability of an output configuration m 
in a realistic network to that of the ideal bosonic network. 
The following result is found. 

Theorem 2.- The difference between the postselected 
probability of an output configuration m in a realistic 
bosonic network and that in the ideal bosonic one satisfies 

|pbd)(^) (26) 

p{n) 

where Ps is given by Eq. m- 

Proof of Theorem 2. - The bound (I26|) follows from Eqs. 
(fT8l) - (l20]) and the following inequality 


E izM\Zr) 

m{l )—rh 


< E (^rl^r) 

m(Z )—rh 


m 

pin) 




(27) 

valid for any Hermitian matrix A whose eigenvalues lie in 
the interval [—1,1]. Setting Aff = |s)(s| — noticing 
that its eigenvalues lie in the interval [—1,1], using Eq. 
(071) . (iroi) . and (OiH) we get 


\p^''‘^\m) — p(rh) 


= (7-Ps 


m{l )=m 

N\ 


< {l-ps)^j^p^‘'^\rh), (28) 

p[n) 


i.e., the bound of Eq. (|26)) . Q.E.D. 

For bosons in pure internal states, as in section 111 Cl 
Eq. (IMl) contains two factors: (i) the complementary of 
the probability of the ideal bosonic behavior 1 — ps = 
1 — per(G)/A^! (see Eq. (IM1) 1 and (ii) the upper bound of 
Eq. (l27l) on output probability p[rh). Recently, a bound 
similar to that of Eq. (l26ll was proved for a real positive 
G-matrix of Eq. ([Mjl and then numerically validated 
for the general case of bosons in pure states . In this 
respect. Theorem 2, besides being an extension to general 
mixed inputs, gives a generally tighter bound than that 
of Ref. [3l|, due to p{n) in the denominator. 

Finally, as was previously noted in Ref. [3l|, the prob¬ 
ability Ps can be read off from the experiment. Indeed, ps 
appears in the output probability (postselected for non¬ 
ideal detectors) of a maximally bunched output configu¬ 
ration (see Refs. [13,113|), say = N. In this case 


N 


n(29) 
T(n) “ P(n) 

and we obtain (details in Appendix [B]) 


p{n) 


(30) 


i.e., the same expression as in Ref. [^, extended here to 
arbitrary (mixed) internal states of bosons. 


III. SMALL ERRORS: SCALING OF THE 
DISTINGUISHABILITY ERROR WITH N 


Let us consider a practically relevant case of small er¬ 
rors, which allows for an analytic solution. We assume 
that detectors are nearly ideal JT = 1 — T 1 and that 
bosons, originated from independent sources, have inter¬ 
nal states which are approximately the same pure state, 

i.e., = na=i®P«, Pa = \4>){4)\ - Spa (note 

that Tr(5pQ,) = 0 and {(flSpalflf) > 0). Below we give all 
expressions to the first order in JT and 6pa- 

Expanding the tensor product inside the trace in Eq. 
dH) we get 







(31) 


with 


Ra=Spa+St\(f)(f>\. (32) 

Let us denote by Gi (cr) the set of fixed points (1-cycles) 
of a permutation cr. Then, computing the trace in Eq. 
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(|6|) for the expression in Eq. (ED we obtain 


= 1- ^ Tr(i?„)- ^ {4>\RM 

CkGCi (o'2<7^^) CK^C'i(o'2<7^^) 

N 

= i-Y^{ci,\sm- Y. ( 33 ) 

“=1 afCi(a2<y^^) 

where we have introduced the fidelity of a mixed in¬ 
ternal state Pa as follows Ta = (ip\pa\<P)- Since the trace 
of the third term on the r.h.s. of Eq. (ESI is zero (all a 
are fixed points for a 2 = <Ji), using Eq. dini we obtain 

N 

Pd = 1 - ^(0|(5f|(/)). (34) 

Of—1 

Therefore, from the definition = J/{pdNl) we have 


Ja 


1 

m 


1-- ^ {l-J^a) 

Pd 

Oi^Cl{cr 2 < 7 i 


(35) 


By Eq. (l20l) the difference between postselected output 
probabilities can be cast as 


p^*'^^(m) — p(m) = (1 — p;,) (m) — p^-*-^ (m) , (38) 


where p*-'*"^ (rh) corresponds to a density matrix or¬ 
thogonal to = |s)(s| of the ideal bosonic input. 

Note that such a j7(-*-)-matrix results in a nonclassical 
probability distribution at a network output. In partic¬ 
ular, = sgn(tT 20 'f^)/N! describes completely in¬ 

distinguishable bosons in an entangled (antisymmetric) 
internal state, simulating the output probability distri¬ 
bution of completely indis ting uishable fermions, due to 
the boson-fermion duality E0|- 

Now, let us consider the maximum (taken over all 
M X M-dimensional networks U) of the trace distance 
between the output probability distribution of an input 
resulting in a density matrix and that of the ideal 
bosonic network 


d = maxZ)(p^*'^\ p*-"'"^). 


(39) 


Erom Eqs. dm and (1551) we obtain 


1 - Ps = 


1 


E E (i--^«)- 

cr a^Ci{(T) 


(36) 


Let us introduce the minimal fidelity = 

min(J^i,..., J^Tv)- Then the bound on the trace 
distance Eq. dm can be estimated as follows 


1 -Ps < 


1-T 

PdNl 


E[iV-|Ci(a)|] 


fV - 1 

-(1-J-), (37) 


where \Ci{a)\ is the total number of 1-cycles and we have 
used the well-known fact that X)cr |C'i('^)l = 1 Sl¬ 
it is seen from Eq. (EZI that a realistic device with 
N bosons is at a constant trace distance to the ideal 
bosonic network for the following scaling behavior of the 
minimal fidelity = 1—0(1/A^). The scaling law derived 
here is weaker by vW then that obtained previously by 
a different method in Ref. [l^. In the next section we 
argue that the upper bound of Eq. dm, resulting in the 
above scaling law, is tight. 


IV. LOWER BOUND ON THE TRACE 
DISTANCE AND DISTINGUISHABILITY 
CONJECTURE 


In section III Al we have already mentioned the prob¬ 
lem of uniqueness of -matrix for a given class of output 
probability distributions. Whereas in all of our calcula¬ 
tions above such a j7-matrix was used as just an auxiliary 
concept whose uniqueness property was never used, now 
we have come to a point when the uniqueness property of 
j7-matrix plays an essential role (though in a formulation 
not as general as in section Hi Al) . 


Then, by Eq. (1551) the maximal over all M x M- 
dimensional networks trace distance becomes (with the 
dependence on parameters given explicitly) 

maxD(p(*^),p) = d(jW,V,M)(l-p,(,7,iV)). (40) 

VC/ 

At the next step the physical distinguishability between 
two mutually orthogonal f7-matrices plays an essential 
role. Let us state this property as the following. 

Special distinguishability conjecture. - There is a con¬ 
stant do > 0 bounding the function d{J^^\N, M) of Eq. 
(1551) uniformly away from zero, i.e., d{J'^-^\ N, M) > do 
for all M, N >2 and all (physically realizable) 77*-such 
that = 0. 

In other words, we conjecture that for any N,M > 2 
one can always distinguish bosons in an internal state 
p(“‘) such that = 0 from the ideal bosonic 

case SnP^™*^Sn = (with output probability distri¬ 

bution of Eq. m) by comparing the respective output 
probability distributions over all unitary networks U. 

If the above conjecture is true, then the upper bound 
of Eq. (1T51) is tight, i.e., it captures the correct scaling 
of the trace distance with N (for instance, the scaling 
of the mutual fidelity for small errors derived in section 
uni is then also tight). E.g., if we take N = 2 the lower 
bound min(d(j7^'’~\ 2, M)) is positive , since is pure 
and corresponds to the ideal fermionic case (see also Ref. 
[3^). 

Erom a physical point of view, at least a restricted ver¬ 
sion of the above conjecture (to some subset of all possible 
J'^^^-matrices), relevant for a realistic BS device, must 
be true. Indeed, it is rather clear that if bosons have 
an internal state which is not a symmetric one, but a 
perturbation of a symmetric state, then the correspond¬ 
ing output distribution must have a nonzero maximum 
distance to the ideal bosonic network. 
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V. CONCLUSION 


Appendix A: Details of the proof of Theorem 1 


We have employed a recently developed approach to 
partially indistinguishable bosons in a multi-mode uni¬ 
tary linear network for derivation of the upper and lower 
bounds on the trace distance of the respective output 
probability distribution to that of the ideal bosonic net¬ 
work. An important property was used, that in an aux¬ 
iliary linear space, where the partial distinguishability 
matrix of bosons is defined, the boson counting measure¬ 
ment at a network output is formally represented also as 
a measurement. This fact allowed us to derive the upper 
and lower bounds on the trace distance by application 
of the standard methods. The form of the upper bound 
has a very clear physical interpretation: it is equal to the 
complementary of the probability of a symmetric input 
state of bosons, i.e., a state resulting in the ideal bosonic 
network. 

In the case of a small distinguishability error there 
is a scaling law of the mismatch of internal states of 
bosons, stating that a realistic device with N bosons is 
at a constant trace distance to the ideal bosonic network, 
if the minimal fidelity of internal states of bosons satis¬ 
fies = 1 — 0{1/N). This scaling law surpasses that 
obtained previously by a different method in Ref. (l^ . 

The derived upper bound is tight (i.e., the lower bound 
has the same scaling in the number of particles as the 
upper one) provided that a physically clear distinguisha¬ 
bility conjecture is true. The conjecture generalizes and 
puts in a mathematical form the universally assumed as¬ 
sumption that a perturbation of an ideal input state of 
bosons generally results in a non-ideal bosonic behav¬ 
ior. This touches upon the very important issue of the 
computational complexity of the BS. Indeed, if the dis¬ 
tinguishability conjecture is false, then there are such 
non-ideal input states of only partially indistinguishable 
bosons, which have the same output probability distribu¬ 
tion as the ideal BS with the completely indistinguish¬ 
able bosons (since the trace distance between the two 
distributions would be zero). 

Our main result can be also states as follows. We 
have given a tight experimentally accessible measure of 
“bosonness” of an unitary linear network with only par¬ 
tially indistinguishable bosons. The measure is given by 
the trace of projection of the internal state of bosons on 
the symmetric subspace of internal states, i.e., by the 
probability of having an input state symmetric in the in¬ 
ternal degrees of freedom (which is also proportional to 
the probability of the maximally bunched output, which 
can be estimated from experimental data). 


VI. ACKNOWLEDGEMENTS 


Let us first verify the following property of vectors | Z-^ 

Y, \Zr){Zr\ = V, IP.,.., = ^ E (^1) 

r 

and show that P is a projector. Indeed, by unitarity of 
a network matrix U 


N M 


= - 75,1111 


aifi) ^ 


^cT2(a) 


N 




M(n) 




Moreover, we have 

-p2 ^ ^ 

(/r(n))2 


E E ^(Tl jTTifT *^0"^ ,7r2'Ti — 7^(7,(T'-; 


<^1 


andtr(7’) = ^. 


Appendix B: On the probability ps 

We need to show that Ps appears in the probability 
of a maximally bunched output configuration, say Ib = 
(?,...,/) (i.e., = N). Recalling that (cr|s) = l/y/N^. 

we get (in this case there is just one measurement element 
for ml^^) 


I m ^ 

= (Bi) 

Thanks to Eq. m, we have saturation of the upper 
bound for an output probability of a maximally 
bunched output in the ideal bosonic case, 

m(^) m(^) 

(B2) 

whereas the other probability appearing in Eq. (1381) van¬ 
ishes, = 0. By resolving Eq. (l38t we obtain 

the result 


p,[n) 


(B3) 
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